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ABSTRACT 


A simple two-level numerical model using the quasi-geostrophic 
forecast equations is developed. Equations are linearized and 
friction is introduced in the surface layer. Solutions are obtained 
numerically by using the initial value approach. Two wind profiles, 
U = aus tanh ye and U = U, anche yly are used and these are known 
to be unstable. For each wind profile the growth rate is determined 
as a function of the wave number. Some observed features of easterly 


waves are reproduced in the numerical solutions. 
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1, INTRODUCTION 


Synoptic scale disturbances in the tropical atmosphere have long 
been recognized (Palmer, 1951; Riehl, 1954). In general these dis- 
turbances possess a "cold-core" structure with the colder air and most 
active weather located east of the wave (Yanai and Nitta, 1968b). 
Vertical motions associated with these waves are generally upward east 
of the wave axis and downward to the west (Yanai and Nitta, 1967). The 
source of energy for these waves is not fully known, however, Chaar 
(1963) argues that when condensation is absent the large scale motions 
tend to be quasi-horizontal and quasi-nondivergent. If this argument 
is correct then the barotropic instability that arises as the result 
of horizontal shear in the Intertropical Convergence Zone could be the 
source of energy for the disturbances (Schminke, 1968). 

Kuo (1949) studied the stability characteristics of a barotropic 
zonal current, and showed that the barotropic zonal current is stable 
if the gradient of absolute vorticity has the same sign in the zonal 
current. Charney (1964) discussed the case where horizontal shear 
instability is produced in the equatorial convergence zone when the 
zone is located away from the equator and the converging air masses 
carry their angular momentum with them. 

The purpose of this paper is to devise and test a numerical model 
that is structured realistically, that is, the perturbation has a con- 
figuration and behavior similar to observed easterly waves. The main 
feature desired is to have the associated weather, and therefore the 


vertical motion, upwind of the trough axis. 


2. BASIS OF THE MODEL 


The assumption that the source of energy for the wave mechanism is 
barotropic instability leads to the selection of a basic current that 
is independent of height. Baroclinicity has been omitted but the 
atmosphere remains stratified. Charney (1963) used a scale analysis to 
show that the barotropic vorticity equation governs tropical motions in 
the absence of condensation. In another paper Charney (1969) discusses 
the role of the phase velocity in energy propagation in an easterly 
regime. Following Schminke (1968) the assumption is made that the 
easterly waves can be roughly described by a two-level quasi-geostrophic 
model. It is recognized that the quasi-geostrophic approximation is not 
very accurate close to the equator, however, the zone of dynamic 
instability is assumed to be associated with the Intertropical Conver- 
gence Zone which is some distance from the equator so that the 
approximation is relatively accurate. 

The type of wind profile chosen will play a large role in determining 
the stability of the disturbance. Jacobs and Wiin-Nielsen (1966) inves- 
tigated barotropic instability and found that there are several unstable 
modes. Since a basic zonal current with no change in the vertical is 
used in this model, the wave with no vertical motion would predominate 
(Jacobs and Wiin-Nielsen, 1966). Easterly waves are observed to have 
vertical motion, so a method must be found to introduce it. This is 
sO8onp Hee by introducing friction in the surface boundary layer. 

Thus vertical motion, which depends upon the stratification of the 


atmosphere, is now forced into the model. 


on 


THE FORECAST EQUATIONS 
A simple two-level model is constructed by dividing the entire 


atmosphere into four layers of constant pressure differential, Ap/2 
(fig. 1), numbered 0 to 4 from top to bottom. 
assumed to be zero at the 


Vertical motion is 


Ap 


Se eS ee re, 
owe i ee 
yg -=<——— en Oe os 
Pars. 1. 
surface is 


Two-level model used for prediction. 


top of the atmosphere, while the vertical motion term at the earth's 
uy = 


- pgw,- 
Charney and Eliassen (1949) used the Ekman theory to derive an 
expression for Wy which is 


(3.1) 


1/2 
1 S. 
= = —- 7 la 
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(22) 


where A. is the kinematic eddy viscosity and a, the inflow angle. 


surface geostrophic vorticity can be approximated by 


The 
~ 2 
Sh og 
Following the development of Thompson (1961), but using different 


(3.3) 
notation, begin with the quasi-geostrophic vorticity equation 
ha 

ot 
Notice that when w 


vy “+ | k x VY¥-V (vy) + 8B 


oo cide 7 
0 Ox o Op ; 
O, the above equation becomes the barotropic 


(3.4) 
vorticity.equation, which is sufficient to describe the main instability. 


Apply this equation at levels 1 and 3 giving 


oy a) 

do _2 2 1 2 

eee lk x Yair 8 5 Semen Sige: (3.5) 
oy (wm, - 0.) 

3 2 2 3 mm 
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Smaller layers could be used if it is desired to stop at some intermediate 
point, say the tropopause, rather than including the entire atmosphere. 


Next consider the quasi-geostrophic first law of thermodynamics in 


the form 
Sale xvy-v CS )+Lo=o, (on 
0 
where 
“ -XT/H 2]. (3.7a) 
pg p 


When this equation is applied at level 2 the result is 


5 CY +¥3) A pa, 
a (tauin Yq Duchalisaaaier ae GAR onic 0. (ARR) 
oO 
Define the following quantities 
ue a a ; (3.9) 
which implies that 
¥, = toot ©, (3.11) 
5 = . A 3 Yap" (3.12) 


Here Yep is proportional to the layer thickness, therefore is a measure 
of its mean temperature. Using these definitions add (3.5) and (3.6) 


and divide the result by 2, obtaining 


oY 


e) m Oy, 
Se vy ate lk xX vy -Vivy ) + lik x vVY-V(WY,) = Bo Se = hi 


o 2p = (. 


(3.13) 
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For the second forecast equation subtract (3.6) from (3.5) and eliminate 





Do using (3.8) 
3 Ds Ze ee. 
a en ay: + |k x VY “VV qu) We Me x VY... 
BY, Wy, 
vivy ) +8. ee oe Gal.) 
where D 
2 aa 
Lb = - (Siel4a) 
Apao 


These are the prediction equations for the model. They can be 
linearized by separating the flow into an east-west current that varies 
only in y, and a small departure from this flow. In this event it is 
possible to treat various waves in x independently, hence the fields 


may be defined as follows 


sa = E(y) + A(y,t) cos kx + B(y,t) sin kx, (3.15) 
Ye = C(y,t) cos kx + D(y,t) sin kx, (3.16) 


where k is the x wave number. 

Substitute expressions for He and Yep thee (3.13) "and (3.14); 
separate the various sine and cosine terms, neglecting all products of 
the quantities A through D, and equating coefficients of the cosine kx 


terms gives 


2 3 
g Gt a) on (BEE? 2h). 


ge -x(25- Jao. (Say) 


For the sine terms the result is 


2 3 
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Repeating the procedure for (3.16) the final equations are 
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S (82 - oe? - Du) = « [= C (k 24,2). mac, 
dy Bey 
ve c| +B Ck + 7c = kK) (Ea oy (3.20) 
Oy oy 
where 
f 8 
K = a 
A ‘ (3225 


It is now necessary to derive the computational equation for vertical 


motion starting with 


pe ee: ae (3.22) 
2 Apo L ot m Ox m dy? : 
where 
Li oY 
ae =e (3.259 
O 
1 oY 
Me = 7) (3224) 


substituting equations @3.23)m(3.24), (3.15) )eand el 6) mco (or ee 
collecting terms, and again neglecting all products of the quantities 


A through D results in 


Wo = mee cos Ss = f AG 1D) 
OD , k_ o£ 
Sin kx ae £ oy y) |. (322) 


This is the equation used in the model for the computation of vertical 
motion. The prediction equations and the vertical motion equation are 
now in the form to be used in the model. 
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4. BOUNDARY CONDITIONS AND FINITE DIFFERENCING SCHEME 


The finite differencing scheme used is illustrated below with a 


sample variable N; 


ON _ 1 
Oy fe thas : Be, (471) 
2 
OMNI 7 
- ye ics 2N, +N, 4) (4.2) 
3 
ON 1 
— = —— N. SNe +N. )- 
Se oH? | Car itl ») 
(N, - 2N. 4 > N._9) | (43) 


where i is the grid index and H is the distance between grid points. 

Centered time differences are used for all quantities except those 
involving friction. The frictional terms are computed at time (t - At). 
The first step in all cases is a forward time step. The second order 
equations for time tendencies are solved by the exact method of Richtmyer 
mo57, p. 101). 

Yanai and Nitta (1968a) studied the problem of finite difference 
approximations in solving dynamic instability problems of non-divergent 
barotropic currents. They showed that for both symmetric and anti- 
symmetric zonal currents the exact boundary conditions can be replaced 
by rigid boundary conditions at a distance equal to the half width of 
the shearing wind belt. At any boundary distance larger than this value 
the numerical and theoretical values tended to have a constant difference. 
Therefore rigid boundaries are placed at y = -w/2 and y = w/2, where w is 
the total width over which computations are made. With the exception of 
E (the basic current) all variables have the boundary condition A = B = 
C=D=0 at y = -w/2 and at y = w/2. When second derivatives are 


required on the boundary they are set equal to zero. 
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Fig. 2. Wind profile U = -U tanh yly 


>. WIND PROFILES AND INITIAL CONDITIONS 


Schminke (1968) investigated and discussed the structure of easterly 
waves from a similar model using one wind profile (U = -U,tanh(y/y )) 
and one non-dimensional wave number (1/2y). It is proposed to expand 
this to include another profile (U = Ujsech’y/y.). Additionally the 
behavior of disturbances at various non-dimensional wave numbers under 
both profiles rill be considered. 

Charney (1963) shows by scale analysis that in the absence of 
condensation, flow in the tropics is governed by the barotropic 
vorticity equation. The assumption is made that the source of energy 
for easterly waves must be barotropic instability. Since potential 
energy is assumed not to be an important energy source the basic 
temperature field may be set equal to a constant and is invariant. 
However, small temperature fluctuations are observed with easterly 
waves, so a two level model is needed (Schminke, 1968). 

If the waves are to grow, the zonal wind profile must be baro- 
tropically unstable. Following the suggestion of Charney (1963) that 
the shear layer near the Intertropical Convergence Zone may be formed 
by transport of angular momentum, the unstable profile chosen to 
approximate this is 

U = -Uptanh iss (a) 
illustrated in figure 2. This seems most appropriate during monsoonal 
conditions. The second profile chosen, 

U=JU seein f (5.2) 
O O 
is probably more applicable to mean flow conditions, and is illustrated 


inetioure 3. 


LS) 
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2 
h-y/y 
l ile U = U sec - 
Wind profi - 
Bae. of 
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Gareta (19560) toundmrenat the prozile given in (5.1) is barotropically 
unstable when O0O<k < 7 Betchov and Criminale (1967) have analyti- 
cally, and by use of the eigenvalue approach, investigated this profile 
throughout the range mentioned above and found the maximum instability 
at a non-dimensional wave number of 0.45 (ky | = 0.45). With these 
results in mind this paper is confined to the range 0.25 to 0.60 for 
the non-dimensional wave number when dealing with the hyperbolic 
tangent profile (5.1). Betchov and Criminale (1967) found a rather 
flat zone of maximum instability for the hyperbolic secant squared 
profile (5.2) in the range 0.80 to 1.05 for the non-dimensional wave 
number. Some computations were done in the lower range (0.25 to 0.60), 
but the majority of the experiments with this profile were confined to 
the range 0.70 to 1.15. Results of experiments with both profiles will 
be discussed in a later section. 

With the wind profile independent of height and with no friction oe 
and Ye should behave independently. Jacobs and Wiin-Nielsen (1966) 
showed that the growth rate for the se field is greater than that for 
the Yep field. The thermal field Yep is zero in the initial state but 
it changes as a result of the Ekman friction term. ie is introduced 
into the model as a random disturbance in y. Introduction of surface 


friction leads to the coupling of the two fields. Thus at time t = 0 


A = random disturbance (y), 65,35 
B = random disturbance (y), (5.4) 
C = 0 (365) 
D = 0 (:,6) 
E = Uy, ln(cosh(y/y)), 7) 
or E = Uy tanh yly (538) 


iby, 


6. COMPUTATIONAL PROCEDURE 


A list of the constants used in the computations is given below; 
y = 200 kilometers, 
w = 2000 or 4000 kilometers, 
U_ = 10 meters per second, 
f =5x Tome per second, 
o = 0.8 meters per second squared per centibar squared, 
B = 0.0ger 2.297x ton per meter per second, 
Ap = 50 centibars, 
A = 10 meters per second squared, 
a = 22.5 degrees, 
u = 2.49 x 10 per meter squared, 
H = 25 or 50 kilometers, 


0.5 or 1G; hours. 


At 
When w = 2000 km and H = 50 km a forecast period of 14 days was used, 
in all other cases the forecast period was 28 days. A sample program 
as developed by Schminke (1968) and modified to fit the IBM 360 computer 
is included as Appendix 1 (note: this program is for w = 2000 km and 
H = 25 km). 

The random disturbance introduced into the ae field, (5.3) and 
(5.4), was generated by a random number generation program, then read 
directly into the computer. The same disturbance field was used for 
all er sarinenes with both variations of 8, H, At, and all variations 
of the non-dimensional wave number. After a few days the disturbance 
adjusted to an exponential growth. For convenience, growth rates were 
computed over the last three days of the forecast period. Once the 


field has adjusted it does not matter when computed, rates are always 
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the same. This was done by assuming the amplitudes could be written as 


x, = Sexp (Un/y.) t,, (6.1) 


1 


xX, = Sexp (Un/y,) to, (6.2) 


where n is the growth rate and S is a constant. By forming the ratio 


Xo/X, > the growth rate can be shown to be 


yoy xX 
ne eae (6.3) 
O 2 i 


Ie) 


7. RESULTS 


A series of experiments was conducted under varied conditions. In 
general the total distance over the grid was 2000 km, but one series was 
done over a total distance of 4000 km. Grid mesh was either 50 km or 25 
km with the smaller distance producing better results. When the grid 
mesh was 25 km a time step of 0.5 hours was used, however the majority 
of experiments were run at t = 1.0 hours. All cases were run with 
8 = 0, then with B = 2.29 x hone meters per second. Many combinations 
of friction, wave number and wind profile were run, and selected ones 
are illustrated in figures 4 through 19. 

Growth rates for the disturbances with both profiles and for various 
non-dimensional wave numbers are shown in figures 4 through 11. Betchov 
and Criminale (1967) found a maximum growth rate of approximately 0.19 
for a non-dimensional wave number of 0.45 using the hyperbolic tangent 
profile. With rigid boundaries, no friction and a fairly large grid 
size (H = 25 km) through the maximum shear zone (width approximately 
400 km), the result (fig. 7) was a growth rate of 0.165 at a wave 
number of 0.45. According to Yanai and Nitta (1968a) the best results 
should be obtained with a minimum of 20 grid steps across the zone of 
maximum shear. Therefore the results for the hyperbolic tangent pro- 
file are considered good, and the trend upward toward the Betchov and 
Criminale (1967) results is evident as the grid size is decreased 
(compare eieces 5 and 7). The inclusion of the 8 parameter had no 
significant effect on the growth rates giving only a slight decrease 


ObeO Ol or less, 
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The hyperbolic secant squared profile results were equally 
encouraging. Just as Betchov and Criminale (1967) lead one to expect, 
a relatively broad zone of maximum growth is found reaching from 0.85 
to 0.95 non-dimensional wave numbers. It is difficult to ascertain just 
what maximum rate Betchov and Criminale (1967) found, however, results 
of these experiments show rates that are significantly smaller than the 
rates for the hyperbolic tangent profile. The inevitable conclusion is 
that the hyperbolic secant squared profile is not quite as unstable as 
the hyperbolic tangent profile. As before, decreasing the grid size 
in the maximum shear area led to larger growth rates. Figures 9 and 10 
illustrate this very well. 

In figures 12 through 15 the amplitude and phase relationships of 
+> Yop and w are depicted for the wave number of maximum growth. 
Amplitudes of ie Yr and w show maxima in the region of greatest 
horizontal shear. The influence of the frictional term is clearly 
seen in the relative amplitudes of ty and 13: The three fields show 
close phase relationships and all phases change rapidly through the 
Shear zone. It appears to be a safe assumption that the temperature 
and pressure fields are in close relation for both levels. The omega 
amplitude (fig. 14) shows a distinct minimum in the center of the 
field but sharp maximum peaks on each side in the maximum shear zone. 
Phase relations between t3 and w might almost be termed ideal with the 
sinking motion north of the shear zone about /0 degrees upwind of the 
amplitude ridge, placing the convergence in the area where it is 
normally observed in easterly waves. 

Addition of the 8 parameter made small insignificant changes in the 


amplitude; phase relationships maintained the same basic pattern as 


2a 


without 8. Variation in the non-dimensional wave number over the range 
0.25 to 0.60 left the pattern basically unchanged while producing the 
growth rate changes discussed above. 

Figures 16 through 19 illustrate the behavior of the hyperbolic 
secant squared profile. Note that 4 far exceeds the other components 
in amplitude (fig. 16) and that the maximum is symmetrical around the 
center of the maximum shear zone. ¥3 is very small indicating that 


friction is probably more effective with the hyperbolic secant squared 


profile. The lower field, ¥,, tends to lead the upper field in west- 


3? 
ward movement. This appears reasonable and is observed in tropical 
disturbances of this nature (Riehl, 1954, pp. 215). This gives the 
entire system a tilt toward the east with increasing height; this tilt 
is not noticeable with the hyperbolic tangent profile. 

The omega amplitude is very sharp peaked and symmetric about the 
Center Maeve VOGT ULES .e aomeaaom is a Yo to the boundaries, are very small. 
In this area the convergence appears to be about 100 degrees out of 
phase with the trough. Near _— 0, along the Intertropical Conver- 
gence Zone, the convergence actually leads the trough. 


Use of a finite value for the 8 parameter made no significant 


difference in the results with this profile. 


Ze 
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Fig. 4. Growth rate for hyperbolic tangent profile, 
w = 2000 kn, H = 50 km, with friction. 
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Fig. 5. Growth rate for hyperbolic tangent profile, 
w = Z2000sim, H = 50 km, no friction, 
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Fig. 6. Growth rate for hyperbolic tangent profile, 
w = 2000 km, H = 25 km, with friction 
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Growth rate for hyperbolic tangent profile, 
w = 2000 km, H = 25 km, with no friction. 
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Fig. 8. Growth rate for hyperbolic tangent profile, 
w = 4000 km, H = 50 km, with friction. 
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Fig. 9. Growth rate for hyperbolic secant squared profile, 
w = 2000 km, H = 50 km, with friction. 
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Fig. 10. Growth rate for hyperbolic secant 
squared profile, w = 2000 km, H = 25 km, 
with friction. 


GSK, 


Zor 
0.80 0730 1.00 EeL0 
KY, 
Fig. 11. Growth rate for hyperbolic secant 
squared profile, w = 4000 km, 
H # 25 km, with friction. 
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Fig. 12. Amplitudes of Y » and ¥ for hyperbolic tangent 


rT ‘1 
profile plotted against arbitrary scale, w = 4000 kn, 
H = 50 km, and wave number # 0.45. 
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Fig. 13. 
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Phase relationships corresponding to amplitudes of fig. 12. 


31 


-3y 
O 


20 100 180 260 


Fig. 14. Amplitude of omega, arbitrary scale for same case as figs. 
l2 and 13. 
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Fig. 15. Phase relationship between Y4 and w, same conditions 
as figs. 12, 13, and 14. 
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Fig. 16. Amplitude of Ya ¥ » and Yr with hyperbolic 
secant squared profile, w = 2000 km, H = 25 kn, 


ky. = 0.95, and scale arbitrary. 
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Amplitude of omega, conditions as in fig. 16. 
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Fig. 19. Phase relationship between omega and ¥3 for conditions 
in fig. 16. 
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8. CONCLUSIONS 


The easterly wave model developed appears to be reasonable if the 
assumption of barotropic instability as the energy source is accepted. 
Basic zonal flow has been introduced by two profiles, both independent 
of height. The temperature variations arose from the action of the 
Ekman friction terms. Waves developed by the model behaved in many 
respects as easterly waves have been observed to act. Results show 
that the growth rate of disturbances is a function of wave number, there- 
fore wave length, as well as other factors. Separation of the 
boundaries to greater distances made some small changes but the 
greater effect is obtained by decreasing the grid size, which also 
required changing the time step size, through the shear Zone. 

Vertical motions produced by the hyperbolic tangent profile are 
very similar to those observed in easterly waves, in that the maximum 
sinking motion is upwind of the lower level ridge and the maximum 
rising motion is upwind of the lower level trough. The disturbance 
maintains its phase relationship in the vertical so that no tilt is 
observed. With the hyperbolic secant squared profile the amplitude 
of the vertical motion outside the maximum shear zone is small and 
again the rising motion is upwind of the trough and sinking upwind 
of the ridge. Within the area of large horizontal shear the con- 
vergent area moves to a position downwind of the low level trough. 

It is not known whether this is observationally verified along the 

Intertropical Convergence Zone. Using this profile, the hyperbolic 
Secant squared, the model has developed a solution that indicates a 
definite eastward tilt with height. This tilt has been verified by 
observations. Perhaps some combination of the two profiles would be 


more realistic in depicting actual waves. 
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A more complete investigation of the effect of friction on the 
symmetric secant squared profile appears warranted. There is, also, 
a definite need to include the effects of condensation in the model. 

Non-linear effects, their role in changing the mean wind profile 
and the structure of the disturbance, are very important and must be 
considered. Use of the primitive equations or the balance equations 
would seem to offer the most effective means of incorporating these 
effects. They appear to be crucial if an easterly wave is to grow 
into a tropical storm. 

Always of utmost importance is the need for more detailed and 
extenSive observations, but nowhere is this as true as of the tropics. 
No model can be developed and adequately tested without detailed 


observations. 


ay 
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